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We examine the influence of remote bands on the tendency toward exciton condensation in a 
system consisting of two parallel graphene layers with negligible interlayer tunneling. We find that 
the remote bands can play a crucial supporting role, especially at low carrier densities, and comment 
on some challenges that arise in attempting quantitative estimates of condensation temperatures. 
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I. INTRODUCTION 

In the last few years, there has been enormous ex- 
perimental and theoretical interest in the properties of 
graphene, a hexaaonally ordered two-dimensional sheet 
of carbon atoms.a Graphene is a gapless semiconduc- 
tor whose valence and conduction bands cross linearly 
near two inequivalent Dirac points (valleys) located at 
the honeycomb lattice Brillouin zone corners. Since the 
Fermiievel in graphene can be shifted by the electric field 
effect, El the density and polarity of the charge carriers can 
be tuned by application of a gate voltage. 

A symmetric voltage applied between the layers of 
a double-layer graphene system can induce an electron 
density in one layer, and an equal hole density in the. 
other. Under these circumstances electrons will tendtrEl 
to form a broken symmetry state in which coherence is 
established spontaneously between separate layers in the 
absence of interlayer tunneling. The broken symmetry 
state is driven by Coulomb interactions between layers 
and favored by nesting between the nearly circular Fermi 
surfaces. In this state the one-particle electron density 
matrix has nonzero components that are off-diagonal in 
the layer indices. A bilayer spontaneous coherence state 
is a type of exciton condensate (see below^i-jwhich has 
particularly interesting transport propertiesuu when the 
two layers can be contacted independently. In graphene 
the possibility of independently contacting two graphene 
sheets separated by a distance Mi several nanometers 
has recently been demonstrated.latl After performing a 
particle-hole transformation on the hole-like layer, the 
broken-symmetry state can be viewed as a spatially indi- 
rect electron-hole pair (exciton) condensate. This is the 
viewpoint we will use throughout this paper. 

In addition to having spin and valley degrees of free- 
dom, electron states in single layer graphene are two- 
component spinors with a sublattice degree of freedom. 
As a consequence, the order parameter for exciton con- 
densation has a four-component sublattice structure. 
The structure of the double-layer graphene system is il- 
lustrated in Fig. |l|, where the blow-up at right emphasizes 
the 4-component structure of the order parameter. In 
the broken-symmetry state all elements of the expecta- 
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FIG. 1: (Color online) Double-layer graphene consists of two 
parallel graphene sheets. The black dots are atoms in sublat- 
tice A and the red dots in sublattice B. The four component 
structure of the condensate is illustrated at right. The exact 
definition of the order parameter A is given in Sec. 
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tion value {c\^Cb^a'), where c\ ^ is the creation operator 
of an electron in the top layer in sublattice a = A,B, 
and Cb^a' is the annihilation operator of an electron in 
the bottom layer in sublattice a' — A, B, can be nonzero 
at each crystal momentum. 



This system was first studied by Lozovik and Sokoliktl 
for kpd ^ 1, with kp the Fermi momentum and d the 
interlayeiji-distajice, and later revisited by a number of 
authors .qEHIHtIj Min et alu evaluated the ground state 
superfluid density of the condensate, demonstrating that 
it is proportional to the carrier Fermi energy, and on 
this basis proposed that the Kosterlitz-Thouless (KT) 
transition temperature for exciton condensation could be 
of the order of room temperature in systems with high 
carrier densities and small layer separations— This sug- 
gestion was countered by Kharitonov et al.fj who esti- 
mated that Tc < 1 mK because the interlayer interaction 
is strongly screened in the normal state. In this paper 
we consider two graphene layers without any dielectric so 
that the effective dielectric constant for interlayer inter- 
actions e = 1. Our work differs from earlier work in that 



we include the influence of remote bands, which play a 
supporting role. When a static screening approximation 
is employed we find that the maximum critical tempera- 
ture is of the order of Kelvins and occurs at rather low 
carrier densities. Since static screening underestimates 
the interlayer interaction, particularly for remote band 
effects, higher condensation temperatures appear to be a 
possibility. 

This paper is organized as follows. In Sec. |l|we present 
the main results of our analysis for the condensate struc- 
ture and the transition temperature Tc. In Sec. |III 
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present our discussion and conclusions. Sec. |V| describes 
important technical details of the calculations used to 
obtain the results described in Sec. ||. 



II. RESULTS 

In this paper we estimate the transition temperature 
and determine the optimal structure for exciton conden- 
sation in double-layer graphene. The details of this cal- 
culation are given in Sec. IV. We consider only order 
parameters that have zero center-of-mass momentum. 
Using mean-field theory we derive a self-consistent gap 
equation, which we solve numerically in two approxima- 
tions. First, to gain physical insight, we approximate the 
screened Coulomb interaction by a contact interaction. 
We determine the strength of the contact interaction by 
performing an angular average over incoming and out- 
going momenta on the Fermi surface. In this version of 
our calculation, we use the full graphene dispersion, but 
we find that substituting the Dirac dispersion with an 
appropriate ultra-violet cutoff gives nearly identical re- 
sults. In an effort to obtain more quantitatively reliable 
results, we have also used a separable approximation to 
the screened Coulomb jjateraction, similar to the form 
used by Lozovik et aZ.Jlil to determine the gap at zero 
temperature. This calculation is carried out within the 
Dirac approximation. The greatest source of uncertainty 
in our work is the approximation used for the screened in- 
terlayer interaction. In all of our explicit calculations we 
use a static approximation which overestimates screen- 
ing. The critical temperatures we obtain will therefore 
tend to be underestimated. 

In the present section we present the results of our 
calculations. We describe the optimal condensate struc- 
ture, which turns out to be the same in both approxi- 
mations. We then discuss the phase diagrams we obtain. 
The contact-interaction approximation leads to unphys- 
ically high values for the transition temperature, as may 
be expected for such a crude model. In contrast, using 
the separable approximation for the screened Coulomb 
interaction, we obtain a Tc of order of Kelvins for typical 
carrier Fermi energies Vg ~ 0.25 — 0.5 mcV and interlayer 
distances d < 4 nm, respectively. 




FIG. 2: Dirac approximation band dispersions for the top and 
bottom layer. The shaded areas correspond to filled states. 
The four components of the order parameter in the band rep- 
resentation are indicated. In the close-band approximation 
only the (lower energy) valence band in the top layer (bottom 
left) and the (higher energy) conduction band of the bottom 
layer (top right) are retained. 



A. Condensate structure 

In the condensed state, the order parameter sublat- 
tice structure (AQ_Q'(k) with a, a' — A, B) is opti- 
mized. (The explicit definition of the order parameter 
A is given in Sec. [V.) We find that the structure that is 
optimal has the form A = [Aaa, ^ab, '^ba, ^bb) — 
(Ai, A2,-A2,-Ai), where Ai > A2. The property 
that Aaa — —Abb can be understood in terms of the 
primary mechanism by which condensate formation low- 
ers the energy of the system, namely the opening of an 
avoided crossing gap between the conduction band of 
one layer and the valence band of the other layer. In 
the Dirac-band approximation the conduction and va- 
lence band sublattice spinors are (l,exp(i0k))/-\/2 and 
(1, — exp(j0k))v2 respectively, where k is momentum 
measured from the Brillouin-zone corner and 4>k is the 
angular orientation of this momentum. An order param- 
eter with Aaa — ~Abb couples these spinors with equal 
strength at all (pk- 

We note that the same order parameter sublattice 
structure is associaied with the broken inversion symme- 
try often thoughttj to be plausible in isolated graphene 
sheets. In that case, of course, the mean-field potential 
couples sites in the same layer. Broken inversion symme- 
try in an isolated graphene sheet is analogous to chiral 
symmetry breaking in elementary particle physics and is 
most likely to occur in neutral sheets without any carri- 
ers. Experiments appear to establish that spontaneous 
gaps do not in fact occur in single-lcw«r graphene; angle- 
resolved photoemission experimentsEj are perhaps most 
conclusive ip_this respect. Spontaneous gaps do how- 
ever appeaitZI in neutral graphene sheets when a mag- 
netic field is applied. It is sometimes arguedllj that the 
appearance of gaps in a field demonstrates that this or- 
der is barely avoided and latent even in the absence of a 
field. (For a contrary view see Ref. [ ll^). We will show 
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FIG. 3: Tc(U) for the contact-interaction approximation for 
Vg =0.1 eV. Here, U is the effective interaction strength in 
eV units. The critical temperature Tc is expressed in units of 
the hopping parameter t = 2.8 eV. The solid line is obtained 
using the full dispersion, the dashed line using the Dirac ap- 
proximation, and the dotted line using the close-band approx- 
imation. 



that, because of their order parameter compatibility, la- 
tent sublattice-pseudospin chiral symmetry breaking or- 
der in a single layer is favorable for bilayer excitonic or- 
der. Although the presence of carriers always acts against 
order in an isolated graphene layer, we show that in the 
bilayer case it acts in favor of order provided that the 
carrier densities have opposite signs in opposite layers. 



B. Full dispersion and close- band approximation 

In this section we consider the influence of the full 
graphene dispersion and the close-band approximation 
on our result for the transition temperature. Since we 
are interested in qualitative effects, it is sufficient to con- 
sider the contact-interaction approximation results. We 
show in Fig. ^ the critical temperature Tc as a function 
of the effective interaction strength U for a fixed car- 
rier Fermi energy V^ = 0.1 eV. The transition temper- 
ature Tc is scaled by the hopping strength t = 2.8 eV. 
Note that the interaction strength U in our model has the 
units of energy and that it is obtained as the continuum 
model interaction strength divided by the unit cell area. 
The solid line is obtained using the full graphene disper- 
sion and the dashed line using the Dirac approximation. 
When we employ the close-band approximation, we ob- 
tain the dotted curve. Below we relate the parameter U 
to the contact-interaction approximation to the screened 
Coulomb interaction: ^'^''(q) — >■ C/, where V^^'^{q) is the 
screened interlayer interaction matrix element. 

When we want to employ the Dirac approximation, we 
are confronted with the fact that the theory has one free 
parameter, namely the value of the cutoff energy ^ that 
is needed to cure the ultra-violet divergence in the gap 
equation on which Tc will depend. In Sec. IV we explain 



how to choose ^, so that the results obtained using the 
Dirac dispersion and the full graphene dispersion over- 
lap for temperatures corresponding to the low energies 
where the two dispersions do not differ noticeably. With 
this choice of i^, we obtain the dashed curve in Fig. | for 
the Dirac model and the solid curve for the full 7r-band 
model. Note that the curves differ only slightly and then 
only for temperatures fc^T > 0.2t at which the difference 
between the full graphene dispersion and Dirac dispersion 
becomes noticeable. Since the results obtained using the 
Dirac and full dispersions do not differ appreciably, we 
use the Dirac approximation for the contact-interaction 
phase diagrams shown in the next section. The dotted 
curve in Fig. B^ is the result of the close-band approxima- 
tion, discussed further below. 

Notice that the critical temperature is weakly depen- 
dent on interaction strength at large U and strongly- 
dependent on interaction strength at small U . The weak 
dependence occurs in the range of values where U is 
strong enough to produce order even when the carrier 
Fermi energy is set to zero, and the strong temperature- 
dependence in the range of U where order is assisted 
by the gate-driven Fermi surface nesting. To understand 
this observation consider the linearized Tc equation which 
is derived in detail in Sec. IV and represented there by 
Eq. dil): 
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Here, v{e) is the isolated layer density-of-states per spin 
and valley per unit cell, nf(e) is the Fermi distribution 
function, and s = — and s = + terms correspond re- 
spectively to close and remote band contributions. (The 
cross contributions between these two pairs vanish iden- 
tically in the linearized gap equation.) When the car- 
rier Fermi energy Vg ^ the two types of particle-hole 
pairs make identical contributions to the right-hand-side 
that decrease with temperature on the scale of the band 
width - 3t and approach l/2t for T ^ 0. (The 1/t de- 
pendence for r — >■ can be understood by noting that 
j^(e) ex e/t^.) The linearized gap equation for Vg — >■ 
has a solution only if U/t > 2. For Vg > the situa- 
tion changes. The remote band contribution decreases 
slightly as these bands are separated from the Fermi 
energy. At the same time the close band contribution 
from small e is enhanced, and in fact diverges logarith- 
mically as T — T' because the density of states v{Vg) is 
finite when the energy denominator vanishes. This is the 
Fermi surface nesting effect. Because the right-hand-side 
diverges, the linearized gap equation will always have 
a solution. On the other hand because the divergence 
is only logarithmic, the critical temperature will be ex- 
tremely small if U is well below the strength at which 
coherence already appears for Vg — ?> 0. This observation 
accounts for the rapid drop in Tc at a particular inter- 
action strength. Because the high-energy contributions 
of close and remote bands are nearly identical, the inter- 
action strength required for a high critical temperature 



is badly underestimated if the remote bands are ignored. 
Note that the lowest temperature illustrated in Fig. y 
- O.OOat is - lOOK. 



C. Phase diagrams 

In this section we discuss the phase diagrams con- 
structed from contact-interaction and separable potential 
approximations. We start with the former to gain physi- 
cal insight and then compare with the latter. In Fig. |[ we 
show how the contact-interaction transition temperature 
versus effective interaction strength curve depends on 
carrier density. From right to left the three curves corre- 
spond to the carrier Fermi energies Vg — (0, 0.2, 0.3) eV. 
In Fig. Q we sees that Tc is an increasing function of both 
U and, for a fixed value of U , also an increasing function 
oiVg. 

This figure provides a nice illustration of the separate 
roles played by the carriers (the conduction band states 
occupied by electrons and the valence band states oc- 
cupied by holes) and higher energy states in forming 
the instability. The contribution of the carriers to the 
linearized gap equation scales with the density-of-states 
at the Fermi level, and hence with Vg. The higher en- 
ergy contribution depends on the density-of-states far 
away from the Dirac point near the model's cut-off en- 
ergy. For V^ = only the remote band contribution 
is present and the system has a quantum critical point 
at Uqcp = 6.25 eV which is indicated in Fig. ^ by the 
dot. Because of the 1/E weighting factor in the gap 
equation combined with the linear in E density-of-states 
of the Dirac model, this quantum critical point satis- 
fies a Stoner-like criterion iy{£,)U = 1, where i/(^) is the 
density of states at the cutoff energy ^. When Vg y^ 0, 
condensation occurs at any value of U, but not at low 
temperatures unless U is close to C/qcp- The carrier 
contribution has a larger relative importance, and val- 
ues of U that are substantially smaller than t/qcp can 
still yield high transition temperatures when Vg is pushed 
toward the largest physically realistic values ^ 0.3 eV. 
The high energy state contribution to the gap equa- 
tion can be captured approximately by replacing the 
interaction among carriers by the effective interaction 
Ucs = C//(l — U/Uqcp)- We note that room temper- 
ature condensation corresponds to kBTc ~ 10~^i, which 
is a low-temperature on the scale of this phase diagram. 



The Tc prediction based on the constant interaction 
model depends on a procedure for constructing an accu- 
rate value of U. We discuss such a procedure in Sec. IV 
where we derive an expression in which U depends on 
Fermi energy Vg and interlayer distance d. We estimate 
U by performing the angular average of the screened 
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FIG. 4: The scaled critical temperature Tc/t versus the ef- 
fective interaction strength U in eV's for several values of 
the carrier Fermi energy Vg for the contact interaction ap- 
proximation. From right to left the curves correspond to 
Vg = (0, 0.2, 0.3) eV. The location of the normal state and 
condensed state are indicated. The quantum critical point is 
marked by a dot. 
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FIG. 5: Tc{d) critical curves calculated in the Dirac approx- 
imation with the transition temperature Tc in Kelvin ver- 
sus the interlayer distance d in nanometer for the contact- 
interaction approximation. From right to left the curves cor- 
respond to Vg — (0.05, 0.1, 0.15) eV. To the right of or above 
a curve the system is in the normal phase, to the left of or 
below a curve the system is in the condensed state. 



interlayer Coulomb interaction over incoming and out- 
going momenta restricted to the Fermi surface. The 
phase diagram thus obtained is shown in Fig. pi where 
we plot the transition temperature Tc in Kelvin versus 
the interlayer distance d in nanometers. From right 
to left the curves correspond to carrier Fermi energies 
Vg = (0.05,0.1,0.15) eV. We observe that Tc is a de- 
creasing function of d, as expected. The behavior of Tc, 
which decreases as a function of V^ , is opposite to the be- 
havior observed in Fig. H in which Vg was varied at fixed 
U. This behavior illustrated in Fig. g occurs because 
screening increases with density-of-states and hence with 
Vg, decreasing the value of U. Both the high transi- 
tion temperatures and the trends illustrated in this fig- 



ure are suspect, however, because static screening at the 
Fermi energy is irrelevant for the high energy states in the 
hnearized gap equation. This conundrum demonstrates 
that rehable Tc estimates are chahenging. 

In Sec. IV we show how to obtain the transition tem- 
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perature Tc as a function of interlayer distance d and 
carrier Fermi energy Vn using a separable-potential ap- 
proximation. In Fig. g we plot the resulting transition 
temperature Tc in Kelvin versus interlayer distance d in 
nanometer for three fixed carrier Fermi energies. The 
solid line corresponds to Vg — 0.1 meV, the dashed to 
Vg = 0.3 meV, and the dotted io Vg ^ I meV. The 
separable potential approximation is more realistic and 
at small d captures the expected increase of Tc with Vg . 
In Fig. we plot the separable potential transition tem- 
perature Tc in Kelvin versus carrier Fermi energy Vg in 
meV for three fixed interlayer distances. The solid line 
corresponds to rf = 2 nm, the dashed to d = 4 nm, and 
the dotted to d = 8 nm. Corresponding points are in- 
dicated by colored markers in the two Figs. || and 0. In 
Figs. we see that Tc has a non-monotonic dependence 
on Vg . This behavior is due to a competition between the 
increase of screening (as in Fig. pi) and the increase of the 
Fermi energy for larger Vg (as in Fig. 0). The maximum 
Tc for fixed d occurs for kpd ~ 0.001 and the height of 
the maximum goes as \/d. This conclusion is, however, 
also based on an approximation scheme that is unreli- 
able for the high energy virtual states that appear in the 
gap equation. We note that in obtaining Figs. pP we 
approximated the polarizibality of Graphene by a con- 
stant, neglecting the increase of 11 with scattering mo- 
menta q > 2fci?- Including this dependence will lead to 
a suppression of the transition temperatures shown in 

Figs. |g. 

As mentioned before, the use of a static screened 
Coulomb interaction is unreliable for high energy inter- 
mediate states in the gap equation. If we assume that our 
static screening estimate U is appropriate for e < 2Vg and 
another estimate U is appropriate for e > 2Vg, the gap 
equation becomes 
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Here Uc is the value of U necessary to solve the Vg ^ 
gap equation: 

where i'{£,) is the density of states evaluated at the ul- 
traviolet cutoff. In Fig. ^ we show the Tc versus d lines 
for several values of U /Uc- From left to right the curves 
correspond to tj /Uc = 0.5, 0.7, 0.9. 
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FIG. 6: (Color online) Transition temperature T^ in Kelvin 
versus interlayer distance d in nanometer for three fixed car- 
rier Fermi energies for the separable-potential approximation. 
The solid line corresponds to Vg — 0.1 meV, the dashed to 
Vg — 0.3 meV, and the dotted to Vt, = 1 meV. The points 
corresponding with Fig. M are indicated by the colored mark- 
ers. 




V„ (meV) 



FIG. 7: (Color online) The critical lines of transition tem- 
perature Tc in Kelvin versus carrier Fermi energy Vg in 
milli-electronvolt for three fixed interlayer distances for the 
separable-potential approximation. The solid line corre- 
sponds to d = 2 nm, the dashed to d = 4 nm, and the dotted 
to Vg = 8 nm. The points corresponding with Fig. H are 
indicated by the colored markers. 



III. DISCUSSION 

In this paper we have discussed a number of possible 
estimates for the critical temperature for exciton con- 
densation in Coulomb coupled graphene layers in which 
a gate voltage has been applied to induce nesting between 
the conduction band Fermi level in the high-density layer 
and the valence band Fermi surface of the low-density 
layer. One of our main results concerns the sublattice 
structure of the inter-layer coherence order, which tends 
to be mainly off-diagonal in sublattice and opposite in 
sign for the AA and BB components. This structure 
is due to the momentum-direction dependence of inter- 
subband phases in both conduction and valence band 
states near the K and K' valleys in graphene. Similar 
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FIG. 8: The critical lines of transition temperature Tc in 
Kelvin versus interlayer distance d in nanometer for Fermi 
energy Vg = 0.1 eV. From left to right the curves correspond 
to U/lJc = 0.5,0.7,0.9. 



observations have been made previouslyJl3E2l These extra 
phases cause contributions to the anomalous inter-layer 
self-energy that are off-diagonal in sublattice to tend to- 
ward small values, yielding approximate cancelation. In 
a bilayer system this self-energy structure leads to com- 
patible self-energy contributions from coherence between 
the bands with nested Fermi surfaces and from coher- 
ence between the two remote bands, the valence band of 
the high-density layer and the conduction band of the 
low-density layer. In addition the sublattice structure 
of bands remains consistent in graphene out to energies 
well beyond the Fermi energies of the nested bands, im- 
plying that contributions to interlayer coherence can be 
expected from high energy states in both carrier and re- 
mote bands. 

Critical temperature estimates for coherent graphene 
bilayers are considerably complicated by the importance 
on the one-hand of screening by the gate-induced carriers 
at low energies and on the other hand of lattice scale cor- 
relations at high energies. For a two-dimensional system, 
mean-field theory overestimates the transition tempera- 
ture to a superfluid and instead a KT-transition should 
be considered. For strong short-range interaction this 
consideration leads to a correction of the transition tem- 
perature by a factor ^ 3.E3 We have reported on a num- 
ber of Tc estimates that are based on momentum and 
frequency independent interactions. When appropriate 
values for the effective interaction strength are estimated 
based on low-energy screening considerations, large crit- 
ical temperatures tend to occur at low carrier densities 
simply because screening is then minimized. These Tc 
estimates likely misrepresent trends as a function of car- 
rier density, since screening at high energies is in fact not 
strongly influences by carriers. On the other hand, they 
do correctly, capture the fact that spontaneous coherence 
would appear even in systems without carriers if inter- 
actions were strong enough. In the contact interaction 
approximation the gap equation in the absence of carri- 
ers implies a Stoner-like criterion in which order appears 



for Uv{£^) > 1, where v{^) ^ l/(ia)^ is the density-of- 
states at the graphene Dirac model ultraviolet cut-off. 
The interaction assumes this form because the density- 
of-states in graphene increases linearly with energy over 
a very broad energy range. It is interesting to compare 
this condition with the corresponding Stoner-like condi- 
tion for density-wave states in a single-layer graphene 
sheet, Uv{(,) > lE3 

In the single-layer case, experiment seems to clearly in- 
dicate that the ground state is not a density-wave state, 
i.e. that U < i/{^)^^. On the other hand, the fact that 
density-wave states appear to occur in the presence of a 
magnetic field may suggest that the criterion is nearly 
met. If so. Fig. || shows that gating-induced Fermi sur- 
face nesting can enhance electron-hole pairing correla- 
tions and induce order at substantial temperatures. 

The difference between our results and the prediction 
by Kharitonov et al.iM3 that Tc is maximally of the or- 
der of milliKelvins deserves further comment. Firstly, 
the value of the dielectric constant for the system we 
consider e — 1 differs from the value e ~4 from Refs. |l^ 
ancO. This difference leads to a drastic difference in the 
prediction for Tc, as mentioned above. Furthermore, in 
their estimate of the upper limit for Tc, Kharitonov et 
al. set the width of the pairing region to zero, while we 
have shown that Tc depends exponentially of the ratio of 
this width to the Fermi energy. This distinction is im- 
portant for the low Fermi energies we consider. A similar 
effect leads to the enhancement of the gap at zero tem- 
perature, as shown in Ref. |1J] . None of these estimates 
include the effect of finite frequencies, which is likely to 
reduce screening effects and raise transition temperatures 
further. 

In conclusion, we have determined the sublattice struc- 
ture of the condensate and discussed the inter-layer co- 
herence phase diagrams predicted by various estimates 
of interlayer interaction strengths. Our calculations 
make no assumptions about the sublattice structure of 
the condensate and take the full graphene dispersion 
into account. When contact interaction approximations 
are used high transition temperatures can occur. For 
the separable-potential approximation, we find transition 
temperatures of the order of Kelvins for the parameter 
range Vg = 0.25 — 0.5 meV and d < 4 nm. This results 
diflers-greatly from the result obtained by Kharitonov et 
al.j^xB that Tc is maximally of the order of milliKelvins, 
mainly because of the important role played in our calcu- 
lations by states at high energy intermediate states with 
e > 2Vg. Careful consideration of the roles of retardation 
effects and coherence in reducing screening could increase 
Tc estimates further. 
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IV. METHODS 

In this section we describe the methods we used to 
obtain the results shown in Sec. pi We use a functional- 
integral approach since it provides a convenient starting 
point to account for non-mean-field effects. 



A. Action of double-layer graphene 

The coherent state path-integral representation for the 
partition function Z is given by 



Z= f d[^P' 



-(Salif>',i>]+Si['4)*,i>])/h 



(1) 



where the integration is over all Grassmann-valued fields 
"0* and f/' that are anti-periodic on the interval [0, ?i/3]. 
The non- interacting action So[■^/J*,^l'] in Eq. (H) is given 
by 

'^ [-^G'o^a,a;a',a'(k,«W„) ?A„.^^/ (k, iw„), (2) 

where the momenta are restricted to the first Brillouin 
zone, the a;„ = 7r(2n-|-l)//i/3 are the fermionic Matsubara 
frequencies with /? = l/ksT the inverse temperature, 
a — A^ B is the sublattice index, and a — t,b the which- 
layer pseudospin index. We ignore the electron spin for 
this moment, to return to it later. In the tight-binding 
model for graphene, the inverse noninteracting Green's 
function from Eq. (kI) is 



Go;a;<T' (k, JW„) 



Sa,a' 



r(k) 



/(k) 



where /i^ is the chemical potential for the a-layer and 
/(k) = -i (1 -H e"*''! -h £-'•'•('"1+""=^)), where t = 2.8 eV 
the nearest-neighbor hopping strength, and ri and r2 
are the lattice vectors E'Ej The interaction contribution 
to the total action Si [^p* , V'] in Eq. (nl) describes the in- 
teraction between electrons via the Coulomb interaction 
and consists of both intra and interlayer terms. The ef- 
fect of the intralayer terms is to renormalize the chemical 
potentials fi^y and horimng strength t and to screen the 
interlayer interactionjljo These terms will be omitted in 
the remainder of this paper. To account for them, we as- 
sume that the renormalized chemical potentials are such 
that in the normal state the top layer is electron and the 



bottom layer hole-like, with equal densities. The corre- 
sponding Fermi levels are denoted by the carrier Fermi 
energy Vg . To incorporate the effect of screening by den- 
sity fluctuations, we use the screened instead of the bare 
Coulomb interaction as the interlayer interaction. Then, 
Si[ip* , ip] consists only of interlayer terms and is given by 

-'0 / / 

r,r' a, a' 

X Ct(i",T)C',6(r',T)^a',b(r',T)V'a,t(r,T), (3) 

where the position summations are over all N unit cell 
positions. Since the transition depends only weakly on 
the precise stacking of the layersjj we make the simplify- 
ing assumption that the stacking is such that the A-sites 
(i3-sites) of the top layer lie directly above the A-sites 
(i3-sites) of the bottom layer. 

The statically screened interlayer Coulomb interaction 
is given by 



T/«"(q) = ^y^"-(i-)e 



Viq)e- 



-qd 



1 - 2V{q)W{q) + (1 - e-29rf) V^{q)\P[q) ' 



(4) 



where 



V{q)^ 



1 27r 



47reoe A q 



is the bare interactionjj eo is the permittivity of the 
vacuum, e is the dielectric constant of the surrounding 
medium which we take to be air with e = 1, d is the 
distance between the graphene layers, and A — 3-\/3a^/2 
the area of the graphene unit cell, with a = 0.142 nm 
the nearest-neighbor distance between the carbon atoms. 
The polarizability 11 is given in the Dirac approximation 
by U{q) = -4i/(Vg) for q < 2kF 

which is the momentum range relevant for the carrier- 
band contribution to the gap equation. The factor 4 
is due to spin and valley degeneracy, ami the density 
of states is v{Vg) = {A/2n)Vg/{3at/2fE3 This form of 
Si[^*,tlj] is valid in the long- wavelength approximation, 
where the interlayer Coulomb interaction is independent 
of the sublattice index. It will be convenient to intro- 
duce the dimensionless momentum variable y ~ qd, so 
that the dimensionless interaction V^'^'^{y) becomes 



v'"iy) 



Vil/d) 



y^"(y/d) 



y + 2j{e)kFd + (1 - e-^v) {-f{e)kFd)^ /y ' 
where we defined 

-qV{q)Ii{q) = l{e)kF = 9.66A:i./e. 



(5) 



Derivation of the effective action 



with the interlayer polarization B given by 



Since the condensed state is a broken symmetry 
state, we cannot resort to perturbation theory to de- 
termine Tc- Instead, we perform a so-cahed Hubbard- 
Stratonovich transformation to obtain an effective ac- 
tion in terms of the order parameter for exciton con- 
densation. Concretely, this procedure entails multiply- 
ing the partition function Z from Eq. (|l|) by a Gaus- 
sian functional integral with value unity over the or- 
der parameter AQ^Q,'(r, r',T) that is on average given 
by V'"'''{y - r')(V'j*„(r, T)^:b,a' (r', t)). By an appropriate 
choice of the parameters in this integral, the interact- 
ing action S'/['0*,'0] from Eq. (0) can be canceled from 
the argument of the exponent in the functional integral 
Eq. P). Then, the integral over the electron fields can be 
performed analytically and an effective action is obtained 
in terms of the order parameter, which is given by 

5'eff[A*,A] = -OTrlog(-G-i) 

X a; „, (k, k + K, zcj,„) A„,,, (k', k' + K, iLu^), (6) 

where oj^ = 2T:m/h/3 are now the bosonic Matsubara 
frequencies, [yksF) (k) is the Fourier transform of the 
inverse interaction in position space l/y'^''(r — r') and 
where G~^ = Gq^ — S with the electron selfenergy given 
by 

ft^ajCTia^cr' (k, ILlI„; k ,UJn') = 

- [(5cr,h'^o-',tAQ'^Q(k',k, iWn - iOJn') 

+6a,tSa'M^*a,a'iK k', iw„' - iW„)] . (7) 

We determine Tc for a second-order phase transition to 
the condensed state by expanding the effective action 
5cff[A*,A] in Eq. (g) to second order in the order pa- 
rameter A* and A. Since we expect the order parameter 
to be translationally invariant in space and time, we only 
consider the contribution of the zero frequency and zero 
center-of-mass momentum components of the order pa- 
rameter A* and A to the effective action S'cff[A*,A] in 
Eq. (ph. For notational simplicity, we define the short- 
hand AQ^Q'(k, k, 0) = AQ^Q'(k). This procedure yields 



k,k'' ai,a2,a3,a4 

^ -''^ai ,Q2 ,CK3 ,CK4 



(k,k')A„3,„,(k'), (8) 






h^/3 



( B, 







s^-^Bi 



;'^Bi B2 \ 



e-"f'Bi Bo e-'^">'B2 -e-">'Bi 
^'^B2 Bo 



-e^-^Si e^'^Ba Bo e'i'Bi 



(10) 



V B2 -e'^Bi e~'^Bi Bo ) 

Here we dropped the k dependence of the Bi, defined 

= arg[/(k)], and 

So = 2B(+,+) + B(-,+) + B(+,-) 

Si=S(-,+) - B(+,-) 

B2 = 2i?(+,+) - i?(-,+) - i?(+,-). 

Moreover 

r,, ^_ ^/(g...«(k))-rv(£.,,fc(k)) 
""^''^''^^ ...(k)-e.,,(k) ' ^11) 

with St, Sfc = ±1 and top and bottom layer dispersions 
s|/(k)|-)-l^g and e,^b(k) = s|/(k)|-y5. 



where 



e«,t(k) 

with s — ±1. The letter B is chosen in Eqs. (10 1^) 
because Eq. (nl|) is the expression of a bubble diagram, 
which describes screening by electron- hole pairs. We note 
that when we perform a particle-hole transformation in 
the (hole-like) top layer, the numerator of the right-hand 
side of Eq. (|ll|) becomes 1 - n/(esj^t(k)) - n/(es2,fc(k)), 
and we obtain the familiar expression for the ladder dia- 
gram from BCS theory. 

Now, we comment on the importance of the real elec- 
tron spin. Each independent fermion species contributes 
to the screening of the Coulomb interaction. Therefore, 
a factor of two due to the spin degeneracy should be in- 
cluded in the expression for the polarizability n((7), as 
we did above. The interlayer Coulomb interaction is to 
a very good approximation independent of spin, and we 
need to consider how the effective action Eq. (||) changes 
if we include the electron spin, in Eq. (0) and Eq. (0). 
The electron spin can be incorporated in our formalism 
by extending the definition of a to include both the sub- 
lattice and spin quantum numbers, so that the order pa- 
rameter has 16 components. The noninteracting Greens 
functions in the expression for B in Eq. (|l^ are diagonal 
in spin, so that the same is true for B in Eq. (^0|) and M 
in Eq. (ffl. Thus, we find that the contributions to the 
effective action of the 4 spin pairing channels are decou- 
pled and thus we find four identical T^ equations, one for 
each channel. Thus, it is correct for the determination 
of the phase diagram to ignore the electron spin in our 
formalism. 



-^^Qi,a2,a3,ct4V 7 ) 
- "T7<Jai,a3 0a2,a4 I T7^ ) (k ^ k )— Ok.k'Bcj^ 



.(k), 



C. Derivation of the linearized gap equation 



The transition temperature T,, is now given by the 
y^> maxinuim temperature for which the matrix M defined 



in Eq. (|[) has a zero eigenvalue, or equivalently the max- 
imal temperature for which for all k, ai, and a2 we have 

k' 03,04 

To get rid of the Fourier transform of the reciprocal in- 
teraction, we multiply with V^" (k" — k) and sum over k 
to obtain the gap equation 



k' ,03,04 

X -t>oi,02,03,04(.'^ j^03,04(.K ). [i-^J 



Below, we use two methods to find approximate solutions 
of Eq. (^2|), namely by modeling the interaction as a 
contact interaction and using the separable-potential 
approximation. 

We remark that the screened interaction in Eq. ( p^ ) 
should in first instance be evaluated at frequency w — 
Ek — ek' • The frequency and wavevector arguments which 
appear in the remote band part of this integral are ones 
for which our static screening approximation is not re- 
liable. It is, however, not immediately clear how to im- 
prove on the approximation we employ because of correc- 
tions to the simple RPA screening function and the role 
of a and a* bands that we do not consider. We there- 
fore choose to use the static RPA screening functions in 
this paper, but remain cognizant of limitations in the 
predictive power of our (or any other) semi-analytic Tc 
calculation. 



D. Approximation 1: Contact interaction 

The contact-interaction approximation is rather crude 
for the Coulomb interaction and can only be used to 
obtain qualitative results for the phase diagram and 
the condensate structure. In this approximation, we 
replace the interaction matrix elements in momentum 
space V^^Ck — k') in Eq. ( |l2[ ) by an effective strength 
U which is an appropriate average of V^^" over its ar- 
guments. It then follows that the components of the 
order parameter Aa-^^a20^) ^^^ independent of momen- 
tum- This approach was used previously for this system 
.t2Hl3 However, because the contact interaction averages 
out the structure of the Coulomb interaction, one may 
only expect to obtain qualitative results using this ap- 
proximation. Setting V^'^'^Ck — k') = U in Eq. (|l^), we 
obtain a 4 x 4 matrix equation 

^01,02 — "77 / , ■t'Oi,02,03,04(,i^ j^03,04 J 



SO that the critical condition is that the 4x4 matrix 
n = {U/N) J2]^ B(k) has eigenvalue 1. We have that 



n 



no,o iiii — iiii 112,0 

u_ I ni,i Ho.o n2,2 -ni,i 

4A^ I -ni,i n2,2 no,o ni,i 

n2,o — iii.i iii^i Ho 



where we defined 

nM = ^E^^(l^)cos{^arg[/(k)]}. 



The eigenvalues of 11 can be computed in closed form. 
Setting the largest eigenvalue — > 1 yields the following Tc 
equation: 



1 = - (2no,o - n2,o - n2,2 



lenf , -t- (n2,o - 02,2) 



(13) 



We remark that the eigenvector corresponding to this 
largest eigenvalue has opposite {A, A) and {B, B) com- 
ponents, a result previously found in a model without in- 
tersublattiee components for the order parameter /S.a b 
and As A.t2l 



1. Dirac approximation 

It is interesting to consider the result for the critical 
temperature Eq. (O) in the Dirac approximation and 
compare its solution to the critical temperature obtained 
using the full dispersion. The linear Dirac spectrum 
is often used as an approximation to the dispersion of 
graphene, where one sets /(k) = {3at/2){kx+iky). Then, 
the Bi (k) depend only on the length of k and it follows 
that the 11, ; vanish for nonzero angular momentum I. 
Then, Eq. (HSh can be written as 



U 



k 

2 Jo „_,-, 



s=±l 



sVa 



(14) 



k', 03,04 



where we used that n/(— e) = 1 — n/(e), a factor 2 was 
added in the second line to account for the presence of 
the valley degeneracy in the momentum integral, and ^ 
is some high-energy cutoff, on which we comment below. 
Again, ^(e) is the density of states for a single valley and 
spin species. The equation Eq. ( |l4| ) describes the situa^ 
tion of band-diagonal pairing, as described previously,llil 
in which there is no pairing between the close-lying con- 
duction bands, and far-laying valance bands. Following 
an approach used to analyze superconductivity in single- 
layer graphene,Ej we can evaluate Eq. (p^ further and 
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obtain 



y r^"" , ,tanh(a;') , , , , x w 

-r+x ax J log(cosn(y — a;) cosn(j/ + a;)) 

^ Jy—x •X 

2x r dx'^-^^^-2\ogcos-h{x), (15) 



with X = /3Vg/2, y = f3^/2 and the dimensionless cou- 
pling constant given by A = Uv{^)/2. The system has 
a quantum critical point only for Vg = and f/qcp = 
6.25 eV where A = 1/2 and C/Qcpi^(^) = 1, which is a 
Stoner criterion for the spontaneous polarization of the 
valence bands in the two layers which are filled for Vn=Q. 
An important point is that the solution of Eq. (n3) de- 
pends on the high- energy cutoff ^, and only when we 
choose a particular value for ^, are we able to compare 
the results obtained using the Dirac approximation and 
the full dispersion. We find this value of ^ by demanding 
that the interaction strength C/qcp at which the quan- 
tum critical point occurs for the Dirac approximation in 
Eq. (15) coincides with the value of C/qcp obtained from 
Eq. (13). This equality leads to the equation 



KO 



'v^i/(k)r 



2N 



which yields S, — 6.83 eV. The results obtained using this 
procedure are discussed in Sec. Q and shown in Fig. H. 
When y 3> X 3> 1 we may approximate Eq. (nSf) to obtain 
a BCS-like result for the transition temperature 



fc,T-:^^IH5exp 



Uv{Vg 



y„ 



c~\ 



(16) 

where C = limj^'_j.oo \ Jq [tanh{x) / x]dx — log(i?) > = 

0.82. Apart from the usual BCS term in the exponent 
— l/Uiy'{Vg), we also find an additional term which scales 
as the length of the pairing region ^ over the Fermi en- 
ergy Vg . This effect will also influence the Tc found in our 
separable-potential approximation and one of the reasons 
that we predict a higher value of Tc for small Fermi en- 
ergies as compared to Refs. |Q and Q. We finally note 
that the gap equation on the close-band approximation 
can be obtained from Eq. ( p^ by only taking the s = — 1 
term of the summation. The close-band results are also 
discussed discussed in Sec. |l| and shown in Fig. 0. We 
note that an equation similar to Eq. ( |lq ) was found in 
Ref. Ijljl for the magnitude of the gap at zero tempera- 
ture. 



2. Estimation of U 

By estimating the effective interaction strength C/ as a 
function of the interlayer distance d, we may transform 
the horizontal axis in Fig. and obtain the phase diagram 
with the transition temperature versus d. We estimate U 



by evaluating the angular average of the screened inter- 
action Eq. (y) over the incoming and outgoing momenta 
restricted to the Fermi surface, k and k', respectively 



Uid) 



L 



ek,ek'=Vg 



dkdk'T/''"(k-k') 



/. 



ek,Ck' 



dkdk' 



= - / d(t)V''"{2kFiiin{(t))), 
which can be easily evaluated numerically. 



(17) 
(18) 



E. Approximation 2: Separable approximation 

In order to obtain a quantitative prediction for the 
mean-field transition temperature, it is not sufficient 
to approximate the screened Coulomb interaction by a 
contact interaction. Instead, we will approximate the 
screened Coulomb interaction by a function which is sep- 
arable in the incoming and outgoing momenta, as was 
also done in Ref. [Q to determine the gap at zero tem- 
perature. We show below how to implement this pro- 
cedure concretely. We will consider s-wave solutions for 
the gap functions of the form 



Aai,a2(2/) = ^auaiivV 



(19) 



where cf) is the azimuthal angle of k, the lax,a2 is the an- 
gular momentum quantum number, and we transformed 
to the dimensionless momenta y = kd. Since we showed 
that the influence of the full dispersion is small, we con- 
tinue in the Dirac approximation, where we may rewrite 
the gap equation Eq. (^ as 

Q:3,Q4 

xB^„„„„3,„,(y')Aa3.o.(y'), (20) 
where we defined the angular averaged interaction 



•^ ai ,02,0:3,04 vy 7 y / / 



2't') 



X e 



' *' V'" ( ^y2 + y,2_2cos(0'-^)yy') , 

(21) 



where we anticipated the fact that the /^ q,/ will be chosen 
such that the right hand side of Eq. (21 ) does not depend 
on (j). We furthermore defined the radial part of B from 
Eq. ( p^ ) as B''^ in the following way 



B 



oi ,02,03,04 



iy/d) = B 



01,02,03 



..Ay)e' 



where B should be considered in the Dirac approxima- 
tion and the values of the integers ".01,02^3.04 can be 
read of from the expression for B in Eq. (|lO|) . Note that 
we did not include an extra factor 2 due to the valley 
degeneracy, since we assume that there is no intervalley 
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scattering. Then, we may argue that each valley pairing 
channel leads to an equivalent gap equation, and we may 
ignore the valley label, similar as we did above for the 
electron spin. The right-hand side of Eq. ( po|) should be 
independent of (jj, which is the case when 



The case / = 



Ia, 



B,B 



Iq, Ia, 



Iq — 1, and / 



B,A 



lo + l. 



for some integer Iq. Then, the gap equation Eq. 
be written in a simplified form as 



^ai.a2{y) 



where 



lie) 



Y. j y'dy'Vi^.^Sv^v') 



vr{y.y')^7r d't^^os^ 



X F^"(Vy2 + j/'2-2cos(0)yy')' 



where V^"{y) was defined above in Eq. (|5|). We defined 
B-'^ as dimensionless form of B^, i.e. with the B{si, S2) 
from Eq. (O) replaced by the dimensionless B defined 
by 



B{st,Sb) = - 



l+oxp[/9'(stj/-fcird)] l+cxp[;3'(si,J/+fcF(i)l 

{sty - kpd) - {.%y + kpd) 



with /3' = (3hvf/d. From these expressions we obtain that 
for a fixed value of kpd the transition temperature goes 
as 1/d and increases linearly with the Fermi momentum 
kp and the carrier Fermi energy Vg. Up to this point, our 
rewriting of the gap equation Eq. (O) is exact, under the 
ansatz Eq. (p^). In order to be able to obtain numerical 
results, we now approximate V^^^(y, y') by a function that 
is separable in y and y' . This approximation amounts to 
choosing a function y^°P such that 



vr{y,y') 



vriyWriy') 



(23) 



From Eq. (g^) we see that the y dependence oi Aa^^^aziy) 
in this case goes as VJ^°p (y) so that it natural to define 

Aai,a.(2/) = Vil'l^^{y)/V^^^^^. The gap equation Eq. @ 
thus becomes a 4 x 4 matrix equation independent of y 



K.,^. = ^ E J y'dy%ZJy')KZJy'^ 



"ai,a2.a3,ai[y )^a3,ai- (24) 
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FIG. 9: We plot Vr{y,y) and [V,"=P(y)]^ by the solid and 
dashed line, respectively, for the cases / = (top graph) and 
I — 1 (bottom graph). 



where y^ is some reference momentum. It is natural to 
choose yl°^ as the position of the maximmn of Vi^^{y, y), 
which yields 



yrJo = and 



yi=i 



2.61-f{e)kFd. 



To gain insight in the effect of this approximation, we 
plot in Fig. I the functions Vriy,y) and [V{"'^iy)f, 
which would fall on top of each other if the approxima- 
tion Eq. (^) were exact. In Fig. ^ we plot V{"{y, y) and 
[yj^''^(y)] by the sohd and dashed line, respectively, for 
the cases I = (top graph) and I = 1 (bottom graph). 

Since [V^''°''(2/)] is always lower than Vf^{y,y) one ex- 
pects that the transition temperatures found in our anal- 
ysis are a lower boundary for the mean-field transition 
temperature for exciton condensation. 



After choosing a functional form of V^^^ we can find 
the transition temperature as the largest temperature for 
which Eq. (BJ) has a solution. For our purposes, it is suf- 
ficient to choose the following form of Vj^'^^ 

vr{y,yf) 



vriy) 



^vr{yT',yT') 
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